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Density functionals for nuclei usually include an effective 3-body interaction that depends on 
a fractional power of the density. Using insights from the many-body theory of the low-density 
two-component Fermi gas, we consider a new, nonlocal, form for the energy functional that is 
consistent with the Fock space representation of interaction operators. In particular, there is a 
unique spatially nonlocal generalization of the contact form of the interaction that preserves the 
p 7//3 density dependence required by the many-body theory. We calculate the ground state energies 
for particles in a harmonic trap using the nonlocal induced 3-body interaction, and compare them to 
numerically accurate Green's Function Monte Carlo calculations. Using no free parameters, we find 
that a nonlocality in the space domain provides a better description of the weak-coupling regime 
than the local-density approximation. 
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Self-consistent mean-field theory is the only practi- 
cal tool to calculate properties of heavy nuclei without 
region-specific parametrizations. Two of the leading im- 
plementations, namely those following Skyrme or Gogny, 
rely on interactions that depend on density, and even 
on fractional powers of density, making them density- 
functional theories rather than Hamiltonian theories. 
We would like to go beyond the density-functional the- 
ory as currently formulated, and possibly back to ef- 
fective Hamiltonian theories, for several reasons. Cor- 
relation energies associated with angular momentum or 
particle number conservation are significant and need 
to be treated outside of the mean-field approximation. 
Restoration of good particle number can be carried out 
in many different ways in a Hamiltonian theory, but 
density-functional theory can lead to inconsistencies if 
the parameterization makes use of nonintegral powers of 
density [l|-[i]. Empirical evidence on nuclear compres- 
sional dynamics favors precisely such non-integral pow- 
ers 0, HI , although there have been attempts to keep only 
integral powers [g, 0|- Also, it is far from clear that fur- 
ther accuracy can be achieved without dropping some of 
the assumptions made, such as the Local Density Ap- 
proximation (LDA). These are our motivations for this 
Letter, to explore nonlocal alternatives to the commonly 
used LDA energy functional. This is not the first time 
the question has been raised; in a 2005 workshop^ the 
problem was formulated: 

"Problem 2. - How can one replace in a nuclear density 
functional: (i) dependence on momentum by dependence 
on density, or (ii) dependence on density by dependence 
on momentum? The fact of life that nuclei arc finite 
systems composed of protons and neutrons must not be 
ignored, forgotten, disregarded, neglected, or otherwise 
assumed irrelevant. The consequences of the proposed 
replacements must be considered in the context of (i) 
constructing functionals from first principles (e.g., how 
to replace the Fermi momentum for the density), (ii) 



conserving symmetries (e.g., how to construct an isospin- 
invariant density functional from microscopic results for 
asymmetric matter), and (iii) restoring broken symme- 
tries." 

Equation of state results for pure infinite neutron mat- 
ter at densities p > 0.04 fm -3 have been commonly 
used to constrain Skyrme and other density-functional 
approaches to large nuclei 0, [Io[ . It has recently become 
possible to use the density dependence of the 1 So gap 
in low-density neutron matter llj to constrain Skyrme- 
Hartree-Fock-Bogoliubov treatments and especially their 
description of neutron- rich nuclei 12j. At low densi- 



ties it is possible to express the ground-state energy as 
an analytically known function of (kpa), the product of 
the Fermi momentum and the s-wave scattering length. 
Thus, finite systems of low density offer a good model for 
testing candidate effective 3-body interactions and com- 
paring with density-functional theory in the local density 
approximation. 

As shown by Lee and Yang [l3j], the ground-state en- 
ergy of a low-density Fermi gas with short-range interac- 
tions can be expanded as a power series in the scattering 
length a. The first 3 terms are: 
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where E/N is the energy per particle and k F is the Fermi 
momentum. 

The corresponding energy density £ expressed as a 
function of ordinary density p (— kp/Sir 2 ) is 
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The first and second terms are just the kinetic and two- 
particle interaction energies of mean-field theory, using 
the scattering-length approximation to the effective in- 
teraction. The third term, which we call an effective 3- 
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body interaction, expresses an energy density functional 
that is proportional to p 7 ^ 3 . 

The origin of that term and its fractional density de- 
pendence may be understood from the graphs in Fig. [T] 
The low-density expansion makes use of the scattering 
length, which is calculated by a ladder sum in the two- 
particle channel. One of the terms in that sum is shown 
in Fig. la. However, in the many-body context, shown in 
Fig. lb, that contribution should be excluded if either £3 
or is below the Fermi momentum. The third term in 
the Lee- Yang expansion is simply subtracting out graphs 
of the form Fig. lb, where the cross on the particle line 
indicates that its momentum is below kp (and thus it is 
not a proper Goldstone many-bod y graph ). In standard 
formulations of many-body theory [14l4l7l | the graph Fig. 
lb is calculated as the integral 
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where n k 
orbital k. 

If we take the second-order perturbative expression for 
the energy, but without assuming Fermi gas wave func- 
tions, and an interaction of contact form the formula is: 
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where <pi r is the orbital function of the i-th particle at 
position r, and occ signifies the fact that the states are 
occupied. 

The minimal generalization of a 3-body contact inter- 
action is one that is a function of two positions, r and 
r'. Such an assumption leads naturally to a precisely 
defined effective 3-body interaction that overcomes the 
problems associated with the p 1 1 3 density functional but 
still reproduces the Lee- Yang functional dependence. 
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FIG. 1: a) A contribution to the effective two-particle interac- 
tion in second-order perturbation, b) The same perturbation 
contribution, in a many-body context with orbitals ki,ka, k$ 
occupied. 



The Fock-space representation of this operator is 



(5) 



The two coordinates r and r' represent the two inter- 
action points in Fig. lb. The prescription to calculate 
expectation values with this operator is to take the con- 
tractions in the Hartree-Fock ground state that corre- 
spond to Fig. lb, where the following triply contracted 
terms are the only ones to survive: 
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Note that since the orbital k 3 is occupied, the graph is 
forbidden (does not occur in a Goldstone expansion) and 
is marked with an x. Its contribution to the effective 
two-body interaction must be subtracted as part of the 
effective 3-particle interaction. 

The function /(r,r') should be translationally and ro- 
tationally invariant, i.e. it can only depend on |r — r'|. 
Moreover, we want the functional that follows from the 
new operator to have the same dependence on the den- 
sity as the last term in Eq. © (i.e. p 7 / 3 ) and we demand 
that it contain h 2 a 2 /m so as to correspond to that term. 
Finally, we see from Eq. ((6|) that /(r, r') should have di- 
mensions of EL 3 . The only way to satisfy all the above 
constraints is with an /(r, r') of the form: 
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To obtain the dimensionless coefficient C, we demand 
that Eqs. ([6]) & ([7]) reduce to the Lee- Yang energy func- 
tional in the limit of uniform matter. In that limit, we 
can express the density matrices as 
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where n k is the Fermi gas occupation factor. Carrying 
out the spatial integrals in Eq. ©, we have 
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We match the result of this integral to the second-order 
Lee- Yang coefficient. This allows us to go back to Eq. 
([5]) and express the 3-body operator as a function of the 
two coordinates r and r' , with an operator structure that 
is guided by the terms in Eq. ©: 
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Since we are dealing with a system of two fermionic 
species, in Eq. (fTO)) we have chosen to call them spin- 
up (t) and spin-down (\) . The two middle terms are free 
to take on either value of the spin. 

A finite range 3-body effective interaction like the one 
we propose in Eq. (|TU)l is much softer than usually as- 
sumed. This fact could have important consequences on 
the density profile of different systems. The price to be 
paid for the conceptual clarity and microscopic deriva- 
tion of our 3-body effective interaction is the emergence 
of challenging computational issues. If N is the number 
of particles under study and N r is the number of ampli- 
tudes in the numerical vector respresenting an orbital, 
then for the usual way of doing Skyrme calculations, 
where the interaction depends only on the one-particle 
diagonal density, the latter is of order NN r . In contrast, 
the form of Eq. (10]) is of order N 3 N 2 . 

With a view to testing the effective interaction of Eq. 
(I10p we consider the model problem of dilute fermions 
in a harmonic trap, as in Ref. 18]. We calculate the 
ground-state energies in density-functional theory with 
the LDA functional, and with Eq. (|10l) replacing the Lee- 
Yang term in the functional. For our purposes, we may 
take the orbitals to be of the harmonic oscillator form 
4>(r) = P[r)e~ vr I 2 where P(r) is a polynomial and v is 
a variational parameter. We write down the variational 
energy of the system using the density as follows: 
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In this expression, the kinetic energy is: 
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and the trap potential energy is: 
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The Si is the second term in the Lee- Yang expansion Eq. 
d2|) and is simply: 
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For the last term, Ejj, we choose two different forms: 
a) one following from the Lee- Yang expansion Eq. ([2]) 
and b) one that corresponds to the new 3-body effective 
interaction given in Eq. (ITU1) . These are: 



Eff = \d\£W = 2h2 f^^ (11 - 2 In 2) / d 3 rpi 



35m 



(16) 



with 

lijk — 



h II - 



d 3 rd 3 r'--^—p % {v,v l )p 3 {v,v')p k {vy) . (18) 



h 2 a 2 



E E E ^ k 

i j k 



(17) 



The equation for E^f contains products of one-body den- 
sity matrices. The sum is taken in such a way that there 
are no repeated terms other than those allowed by the 
restriction (clear by inspecting Eq. (fT0|)) that only the 
tt-l an d til configurations are allowed. 

We examine a system of 8 particles (4 spin-up and 4 
spin-down). We choose 8 particles since this is the small- 
est non-trivial closed shell system: we find 128 terms in 
total (64 for fU and 64 for UD- In Table U we show 
the results of variational minimizations of the function- 
als given in Eq. (|12[) in varying levels of sophistication. 
Shown are the energies of the system at the minima. We 
list the results when we keep a) only the kinetic energy 
and the external potential terms in Eq. (|12|) , b) the same 
two terms plus the next one (£/), which is proportional 
to the scattering length a, c) the above three terms plus 
the Lee- Yang highest order term (£fy y ), proportional to 
a 2 , and d) similarly to the previous case but with E\J 
replaced by Eff . 

To test the accuracy of the different approaches, we 
have performed Green's Function Monte Carlo simula- 
tions, which have already been proved to be dependable 
in describing the electron gas, light nuclei, and cold atoms 
23|. As before, we study 8 trapped fermions, assum- 



ing a Hamiltonian of the form: 



^ = E(-^ + ^ 2 ^) + E^') 



k=l 



(19) 



i<3 



where N is the total number of particles and co is the 
trap frequency. The interaction between the particles is 



TABLE I: Results for the ground-state energy (in units of hw) 
of a harmonically trapped system of 8 pairwise repulsively 
interacting particles. The interaction is at the dilute limit, 



E [hw] 


Kinetic + External Energy 


18.0 


DFT (LY LO) 


19.197 


DFT (LY LO + NLO) 


19.436 


DFT (LY LO + H 3 NLO) 


19.465 


GFMC 


19.485(1) 
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FIG. 2: The ground-state energy (in units of Jko) of a system 
of 8 particles plotted as a function of the scattering length (in 



units of y/h/rruj) 



taken to be of the modified Poschl- Teller type: 
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where fi is a parameter that is related to the inverse of 
the effective range and vq is a parameter we can use to 
tune the s-wave scattering length. The scattering length 
for Table U is a = 0. 2082 y/K/muj. For that case, the nu- 
merical calculations were performed taking vq = 10.28 
and pT 1 =0.1 in units of the oscillator length y/h/muj. 
Since the effective range is much smaller than the oscil- 
lator length, this interaction is appropriate to describe a 
low-density system. Also, we have taken the interaction 
to be repulsive to avoid problems with superfluidity. We 
believe that the errors associated with the fixed-node ap- 
proximation are inconsequential at the level of accuracy 
we are considering here. 

We are now in a position to compare the various ener- 
gies from TableUto the microscopic results. The Green's 
Function Monte Carlo result is close to the density- 
functional theory calculation which includes all the terms 
in the Lee- Yang expansion, but even closer to the re- 
sults following from our new term H3. We have also 
extended these calculations to different couplings, and 
gather our results in Fig. [2] At weak coupling, the new 
operator leads to results that are identical with the Lee- 
Yang a 2 dependence, while as the coupling gets stronger 
the i?3 contributions are more repulsive than in Lee- Yang 
(though they have the same power- law behavior), and 
thus provide a more accurate description of the micro- 
scopic simulation. Importantly, our approach contains no 
free parameters: Eq. (|10[) contains the parameter C, but 
that is matched to the appropriate Lee- Yang coefficient 
for uniform matter, while our subsequent calculation was 



performed for a finite particle number in a harmonic trap. 

In summary, we have attempted to combine an aware- 
ness of the theory of the weakly interacting 2-component 
Fermi gas with the desired behavior of density function- 
als, and have proposed a new form of the effective 3-body 
interaction that has a finite range. By including a new 
term of the form of Eq. (|10[) . we are providing realis- 
tic density functionals with a way to match the analyt- 
ically known behavior at very weak coupling. We want 
to stress that to implement the proposed 3-body opera- 
tor and functional what is required is only the one-body 
density matrix. This is computationally more challenging 
than nuclear physics functionals based on the local den- 
sity, but it is still much simpler than approaches (such as 
coupled-cluster theory or the MBPT approximation of 
quantum chemistry) that deal with correlations explic- 
itly. 
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